Dissipation-boosted entanglement in coupled harmonic oscillators 
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We investigate the effect of dissipation on indirectly coupled harmonic oscillators and show that the 
degree of entanglement of initially separable states, caused by squeezing, is enhanced by dissipation. 
This counterintuitive enhancement is generic, since it arises generally when squeezing is present, 
and squeezing is always present when the rotating-wave approximation is not made. The effect 
vanishes if all the oscillators are coupled to identical heat baths, suggesting that "heat flow" may be 
a necessary condition for dissipation-boosted quantum correlations in coupled harmonic oscillators. 
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One of the strangest features of quantum mechanics, en- 
tanglement, is at the same time a crucial resource for 
quantum information processing [ll-Q- It is wei l known 
that superposition states, including entangled states, are 
extremely sensitive to noise and dissipation. For exam- 
ple, decoherence induced by an environment tends to 
reduce quantum coherent superpositions to incoherent 
mixtures Q, an effect one typically wishes to minimize. 
Somewhat counterintuitively, however, decoherence can 
also be used to generate entanglement Schemes to 

that effect, typically involving atoms coupled to cavity 
fields, can be modeled as few-level quantum systems cou- 
pled to harmonic oscillators. In this Letter, in contrast, 
we are motivated by recent proposals to generate entan- 
gled states of micro- or nanomcchanical oscillators 
We investigate the effect of dissipation on a set of pair- 
wise linearly coupled harmonic oscillators, in a setting 
which does not involve any few-level systems. 

Our aim is to investigate whether dissipation could be 
used to generate or enhance entanglement among cou- 
pled quantum harmonic oscillators initially prepared in 
"classical" separable Gaussian states, such as coherent 
or thermal states. Such a scheme could then be applied 
e.g. to micro- or nanomechanical oscillators coupled to 
cavity fields, and would be analogous to existing schemes 
involving quantum harmonic oscillators coupled to few- 
level systems. One easily realizes, however, that it is 
not possible to generate entanglement from classical ini- 
tial states of quantum harmonic oscillators coupled by 
only passive, i.e., nonsqueczing linear interactions. This 
is because the state of any number of passively coupled 
harmonic oscillator modes that are initially in a classi- 
cal state will remain classical. This result still applies 
in the presence of decoherence and noise since this can 
be modeled via passive coupling to additional harmonic 
oscillator modes. Essentially, what one has, both in the 
presence and absence of dissipation, is equivalent to a 
linear optical network involving only beam splitters and 
phase shifters, which cannot generate any entanglement 



starting from classical states [14[. One therefore con- 
cludes that if all couplings are passive, the presence of 
few-level systems is crucial for dissipation to be able to 
enhance entanglement. 

The above equivalence with linear optics, how- 
ever, relies on making the rotating-wave approximation 
(RWA) [l5(, which can profoundly affect entanglement 
properties. For example, if one makes the RWA then 
the ground state of a chain of harmonic oscillators is 
separable, while if the RWA is not made, it is entan- 
gled [1(1]. More generally, if the RWA is not made for 
linearly coupled oscillators, then generic squeezing inter- 
actions remain and the coupling is no longer passive. In 
this case entanglement between the oscillators, initially 
prepared in "classical" states, may indeed arise. We have 
found that, remarkably, decoherence can then act to in- 
crease the degree of entanglement. This appears to be 
a stochastic resonance-like phenomenon, reminiscent of 
that found for a cavity QED setting Q, and it may be 
a general phenomenon. Interestingly, it seems that en- 
tanglement is enhanced only if the oscillators are coupled 
to baths with different temperatures and/or dissipation 
rates. 

In this Letter we investigate in detail two indirectly 
coupled quantum harmonic oscillators, where the inter- 
action between the two oscillators is mediated via linear 
couplings to a third middle oscillator. We will also re- 
mark briefly on the scenario in which the two oscillators 
are directly coupled to each other. 

Denoting the vibron excitation creation (destruction) 
operators for the two end oscillators by a^(a) and b'(b), 
and for the mediating oscillator by $(c), and assuming 
that all three oscillators are resonant, the Hamiltonian is 
given by 

H = hoj(a^a + tfb + c^c) + hK(a + a^ + b + tf)(ci +c), (1) 

where k is the coupling strength between end and mid- 
dle oscillators. We assume that each oscillator is coupled 
to its own independent equilibrium heat bath. The time 
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evolution of the system in the Born-Markov approxima- 
tion is then described by the Lindblad-typc master equa- 
tion 



Recalling that A(£) is symmetric, we can write 



dp 

dt 



[H, p] + jC a p + -fCbP + T(n + l)C c p - 



TnCjp, 
(2) 

where 7 is the decay rate of the two end oscillators, 
coupled to zero-temperature heat baths, V is the de- 
cay rate of the mediating oscillator, which is coupled 
to a heat bath with average thermal occupancy n, and 
C x p = 2xpx' — x'xp — px)x. 

We shall solve the master equation ([2J using the full 
Hamiltonian ([T]). For later reference we note that the 
RWA amounts to dropping terms with two annihilation 
or two creation operators, so that the interaction part of 
the Hamiltonian ([1]) reduces to 

#rwa = ft«(ac t + ca) + be) +cP). (3) 

The RWA is a good approximation if k <C oj. It greatly 
simplifies analytical and numerical treatment, but fails to 
capture the quantum correlations that arise for classical 
states evolving according to ([1]). 

Analyzing the entanglement dynamics for continuous 
variable systems is rather involved unless one only con- 
siders Gaussian states, which include thermal, squeezed 
and coherent states. This fits our purpose, since we con- 
sider classical initial states, and the Hamiltonian con- 
tains squeezing terms but no higher order terms. To 
solve the master equation with the full Hamilto- 
nian ([T]), we define a normal-ordered characteristic func- 
tion [H] as x(u,v,w) = (e" al * e~ u ' a e v ' b ' * e~ v *' b e w& ' *e~ w ' & ). 
From the characteristic function one can obtain the ex- 
pectation values of quantum mechanical observables, e.g. 

(at™(t)6t»(t)) = (&r(&rx(«, v, w = o, t)\ u , v=0 . _ 

If the oscillators are initially prepared in Gaussian 
states then it is consistent to make a Gaussian ansatz for 
the time-evolved characteristic function, x(it, v, uj, i) = 
exp [— z T A(t) z + iz T b(t)~\ . Here A(£) is a time depen- 
dent 4x4 symmetric matrix, b(t) is a 4x1 time de- 
pendent vector and z T = (x\, X2, W2, Wi) with xi = 
(u + u* + v + t>*)/4, X2 = (u — u* + v — v*)/Ai and 
w = wi + iw2- The corresponding partial differential 
equation for x( u i v i w i t) then becomes [IH 
d 

■ 2kz t Nz)x(u, v, w, t) = /MVx(n,!),ii),i), (4) 



where V = ( 
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(5) 
(6) 



Here is the 2x2 null matrix, I is a 2 x 2 identity matrix, 
and 

' ° M v=( 4k °) V=(° ° 
-10 ' loo/' I 2k 



s = 



(7) 



A(f) 



P(t) Q(t)\ 
Q(t) T R(*) J 



(8) 



where P(t) and R(i) are 2x2 symmetric matrices. Using 
the above form of A(i), Eq. ((4]) reduces to the three 
coupled matrix differential equations 



P 
R 

Q 



w(SP - PS) - 2 7 P + UQ T + QU (9) 
w(RS - SR) - 2rR + VQ + Q T V + 2nn 
w(SQ + QS) - (r + 7)Q + PV + UR + 2kI, 



along with b(t)=M.b(t) . This set of equations can be 
numerically solved for the time dependent matrices P(t), 
R(£) and Q(t), which together with the solution b(t) = 
exp(M£)&(0) allows us to obtain the two-mode quantum 
characteristic function x(it,i>, w = 0, t). For the steady 
state, which is the regime of current interest to us, an 
analytic expression for )((ii,d,ii) — 0, t) can indeed be 
derived but it is rather lengthy and impenetrable. 

Noting that an initial Gaussian state evolving accord- 
ing to Eq.(j2|) will always maintain its Gaussian char- 
acter, allows us to fully characterize the quantum cor- 
relations between the two end oscillators in terms of 
their Wigner covariance matrix. In the steady state, 
the covariance matrix V is a 4 x 4 stationary sym- 
metric matrix with V%j = {{RiRj + RjRi))/2 where 
i,je{a, b} and R T = (q ai p ai q b ,p b ). Here & and pi are 
the position and momentum quadratures of the ith os- 
cillator. A widely used entanglement measure is loga- 
rithmic negativity, which is an entanglement monotone 
and easy to compute. For a two-mode Gaussian con- 
tinuous variable state characterized by the covariance 
matrix V, the logarithmic negativity is computed as 
N = Max[0, -\og{2vJ)\ where 2/_ is the small- 

est of the symplectic eigenvalues of the covariance ma- 
trix given by v_ = ^Ja/2 - ^[{o 2 - 4DctVy/2. Here 

Ai Cj 



a = DctAi + DctBi - 2DetCi and V 



Br 



where Ai (Bi) accounts for the local variances of mode 
a(b), and Ci for the inter-mode correlations. 

The steady state value of the logarithmic negativity for 
the two coupled oscillators was evaluated and is shown 
in Fig. [T] Remarkably, steady state entanglement be- 
tween the two coupled oscillators is a non-monotonic 
function of the dissipation strengths T and 7. For a 
given dissipation strength 7, there is an optimal value 
of r for which entanglement between the two indirectly 
coupled oscillators is maximized. This counterintuitive 
effect is the main result of the paper and is a clear sig- 
nature of the possibility of increasing the quantum cor- 
relations by purely incoherent means. Some insight into 
why a certain level of dissipation may enhance the en- 
tanglement can be obtained by transforming to the in- 
teraction picture, where the time evolution of operators 
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FIG. 1: (Color online) Logarithmic negativity Af, quantify- 
ing the entanglement between the two indirectly coupled har- 
monic oscillators in the steady state, plotted as a function of 
r and average thermal occupancy of the bath n where u = 1, 
ft = uj/3 and 7 = cj/10. 



Oi(t) = cxp(iHt/h)Oexp(—iHt/H) is governed by the 
Hamiltonian ([1]). In this way the entangling effect of 
the Hamiltonian part of the evolution is separated from 
that of the dissipation. Without dissipation, the state 
pi does not evolve and the entanglement induced by the 
Hamiltonian part of the evolution is manifested in the 
time-dependence of &i(t), bi(t) and ci(t). In the interac- 
tion picture, including dissipation, the state pi(t) evolves 
according to 

j tPl (t) = [ 1 £ I a + 1 C I b +T{n+l)C I c +Tn£ I c ,]p I {t), (10) 

where C a ,C^,£, c and C 1 ^ are obtained by replacing the 
time independent a,b,c,c f with ai(t),bi(t),ci(t),c\(t). 

A full diagonalisation of the Hamiltonian in ([I]) will 
be given elsewhere (H)]. The essential feature is that 
creation and annihilation operators are coupled during 
the time evolution, that is, 

[a T (t) , h (t) , ci (t) , a\ (t) , b\ (t) ,c\(t)] T = 
F.[a(0),&(0),c(0),at(0),6t(0),ct(0)] T , (11) 

where F is a timc-dependent 6x6 matrix. As a result, 
the dissipation terms in the interaction picture will con- 
tain not just terms with one creation and one annihilation 
operator, but also terms with two creation or two annihi- 
lation operators. This may give rise to additional squeez- 
ing. The noise may therefore generate further entangle- 
ment over and above that generated by the Hamiltonian 
H in (flj. This explanation also indicates that enhance- 
ment of entanglement through dissipation should be a 
general phenomenon for non-passive linear interactions, 
since whenever squeezing is present in a Hamiltonian, the 
dissipation terms in a similar interaction picture will in 
general contain additional squeezing terms. 

On the other hand, no boosting of entanglement is 
expected for isotropic dissipation. This is the case for 
instance when all oscillators are coupled with identical 
decay rates to baths with the same nonzero temperature 
n. The no boosting argument relies on the property 



where s,f,g refer to the independently evolving nor- 
mal modes for the Hamiltonian ([1]) [l9[, with s — (a — 
b)/y/2, f = {a + b)/2 + c/y/2 and g=(a + b)/2- c/y/2. 
Thus, isotropic dissipation for the a, b, c modes corre- 
sponds to isotropic dissipation for the independently 
evolving normal modes s, /, g. Since such dissipation can 
only degrade any squeezing of the normal modes, where 
this squeezing is expected to correspond to enhancement 
of entanglement between the a and b modes, asymmetry 
between the reservoirs therefore seems necessary for dis- 
sipation to enhance quantum correlations between har- 
monic oscillators. This is intriguing, since it indicates 
a possible connection between enhancement of quantum 
correlations and heat flow. 

When the RWA is made, the dissipation terms in 
the corresponding interaction picture only contain terms 
with one creation and one annihilation operator, in con- 
trast to Eq. (fTT|) . and thus are not expected to create 
squeezing. This is consistent with the general argument 
that no entanglement can be generated if all couplings 
among system and environment modes are passive. An 
easy way to see this is to consider the Glauber-Sudarshan 
P-representation 15|. Any two- mode state P12 can be 



written in terms of a function P{a, (3) as 



P12 



d 2 ad 2 /3P(a,/3)|a)n(a| 



122 



09|. (13) 



For a classical state, the P function is positive, and 
the state is separable. If the P function is negative, 
however, this is a signature of non-classical correla- 



tions in the system |21[. Coherent states are always 
transformed into other coherent states when the cou- 
pling is passive and linear [l5| . For example, using a 
50/50 beam splitter, \a)i <E> |/3)2 may be transformed into 
I (a + /?)/\/2)i (8 I (a - /3)/V2) 2 . It therefore follows that 
under such couplings an initially positive P function will 
remain positive, meaning that a classical state will re- 
main classical and separable |14J. This is true also when 
dissipation is included, since dissipation can be modeled 
via passive linear interactions with a reservoir [22| . 

Within the RWA an exact analytical solution for the 
master equation is also possible [19j . For example, if the 
two end oscillators with 7 = and n = are initially 
prepared in thermal states with mean occupancy n', the 
steady state is a thermal mixture of "dark states" , 



P = 



E ^ + Dn+l l^nX^nl^lOcXOel, 

71—0 ^ ' 



(14) 



(C a + C b + C c )p = (C s + Cf + C g )p, 



(12) 



where |*~) = (s t )"|0) a |0) 6 /(V^I) and s = {a - b)/y/2. 
The dark states |^^) are so called since they are not 
affected by dissipation in mode c. Each dark state indi- 
vidually is entangled, but interestingly, the mixture in 
Eq. (Tnj is nevertheless separable, which can be con- 
firmed by checking that its P function is nonnegativc. 
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We note that the RWA Hamiltonian ([3]) looks very sim- 
ilar to the one studied in Q , where it is found that noise 
does enhance steady-state entanglement even if the RWA 
is made. The main difference is that in Q, two cav- 
ity modes were interacting with a noise-driven two-level 
atom, leading to entanglement between the two cavity 
modes. In our case the atom is replaced by a third oscil- 
lator, and as we have argued above, an initially classical 
state of any number of passively coupled harmonic os- 
cillator modes will always remain separable if the RWA 
is made. The presence of few-level quantum systems 
is therefore necessary if entanglement is to be created 
with only passive linear couplings, whether dissipation is 
present or not. 

It is interesting to note that nonclassical correlations 
beyond entanglement, as quantified by Gaussian quan- 
tum discord [20(, show similar behavior as entanglement. 
Details of this will be provided elsewhere [l9|, but the 
main result shows that these quantum correlations sur- 
vive in the steady state. 

Our results may be relevant, e.g., for a chain of coupled 
nano- or micromechanical oscillators 0] , for two nanocan- 
tilevers indirectly coupled to an ultracold Bose-Einstein 
condensate (Hf , or for a chain of trapped ions in a Paul 
trap [lH • Direct Coulomb interaction between three such 
ions can provide an interaction of the form in Eq. (fTj). 
Our study shows that dissipation can play a construc- 
tive role in enhancing quantum correlations in continu- 
ous variable systems, and suggests a general method for 
checking if this might be the case. We have also stud- 
ied the dynamics of only two directly coupled harmonic 
oscillators, both of which are asymmetrically coupled to 
their independent reservoirs. This physical scenario is 
not exactly the same as discussed before. This is because 
the ground state of a chain of linearly coupled harmonic 
oscillators can only have nearest-neighbor entanglement 
Q. In this case, too, we have confirmed the constructive 
role of local dissipation in maximizing quantum entangle- 
ment of coupled harmonic oscillators in the steady state, 
although the effect is not as unambiguous as with the two 
indirectly coupled oscillators described in this Letter. 

In conclusion, we have studied the dynamics of two 
indirectly coupled harmonic oscillators initially prepared 
in classical Gaussian states, with interactions mediated 
via a third oscillator. If the rotating-wave approxima- 
tion is not made, then quantum correlations are induced 
between the two end oscillators. Strikingly, the degree 
of quantum correlations is maximized for some optimal 
value of the dissipation rate V . The connection between 
heat flow and the enhancement of quantum correlations 
merits further investigation. 
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